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Abstract: The 2-matrix model can be defined in a setting more general than polynomial 
potentials, namely, the semiclassical matrix model. In this case, the potentials are such that 
their derivatives are rational functions, and the integration paths for eigenvalues are arbitrary 
homology classes of paths for which the integral is convergent. This choice includes in particular 

the case where the integration path has fixed endpoints, called hard edges. The hard edges 
induce boundary contributions in the loop equations. The purpose of this article is to give the 

loop equations in that semicassical setting. 



1 Introduction 



X 

Orthogonal polynomials and biorthogonal polynomials, in the context of random matrices, have 
been mostly studied for polynomial potentials, on the real axis, or sometimes on homology classes 
of contours going from oo to oo ^H|- However, it is possible to define matrix models corresponding 
to a more general context, in particular, the "semi-classical" (called so, because it contains all the 
classical polynomials). It is defined as follows (see Bertola [l]): 

Consider two potentials, (i.e. functions of a complex variable), V\(x) and 14(2/) whose deriva- 
tives V[(x) and V^iy) are rational functions (notice that oo maybe a pole of V{ (resp. V£), this is 
the case if V[ (resp. V£) is a polynomial). 

Then consider a generalized integration path Y = Yliij K i,jli x 7i> sucn that the following 
integral is absolutely convergent: 



e 



.K mx)+V2{ y)- X y] dxdy 



r 
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The possible paths 7 (resp. 7) are described in jT], they can be closed or open: 

- If 7 (resp. 7) is a closed contour, the result is non-zero only if it encloses a singularity of 
e - ^ 1 (resp. e~T V2 ). 

- if 7 (resp. 7) is an open contour, its extremities can be: 

• any point in the complex plane, except at the poles of V[ (resp. V 2 ), 

• a simple pole of *jVl (resp. yV^), with residue G R_, 

• a degree > 2 pole of V[ (resp. V 2 ), if 7 (resp. 7) approaches the pole in a sector where 
MfVi > (resp. V 2 > 0). 

Some of the extremities are such that e~^^- Vl ^ +V2 ^~ xy ^ vanishes, they are the poles of V{ 
(resp. V 2 ) of degree at least 2, as well as the simple poles of (resp. yV^'), with residue G M_, 
we will call them "poles". And some extremities are such that e^TW^+^fe) - ^] does not vanish, 
they are arbitrary points in the complex plane, we call them "hard edges" , and write them: 

(X u X 2 ,...,X Kl ) (resp. (Y 1 , Y 2 , . . . , Y K ,)) (1-2) 

Definition 1.1 Consider a path 7 in the complex plane. Define the set of "normal matrices" on 
the path 7: 

H N (<y) := {M G GL N (C) \3U G U(N), 3(x u . . . , x N ) G 7^, M = Udiag(x u . . . , x N )U^} (1-3) 
with a measure: 

dM := [ Y[(xi - Xj) 2 J dU f[ dx { (1-4) 

\i>3 J i=l 

where dU is the Haar measure on U(N), and dxi are the tangent vectors to the path 7. For a 
given M G Hn(T), the matrix U and the eigenvalues Xi are not uniquely defined, they are defined 
up to a permutation, and a conjugation by a diagonal unitary matrix, however dM is invariant 
under these operations, and thus is well defined. 

For example, if 7 is the real axis we recognize if a? (7) = if at the set of hermitean matrices. 
Then, as in [T2|, we define the set of normal matrices on a generalized path T = ^\ . Kijji x 7,-: 

Definition 1.2 Let F = Ylij K i,j1i x lj be a linear combination of products of paths. We define: 

H 2 N (T) := {M h M 2 G GL N (C)2 \ 3£7 ls U 2 G U(N) 2 , 3(x u ...,x N )e C N , 3(y x , ...,y N )e C N , 
Mi = C/idiag(xi, . . . , x N )U\, M 2 = U 2 dia,g(y 1 , . . . , y N )U\, 



Vfc = 1, . . . , N, 3i k x k G 7 ifc , Vfc = 1, . . . , iV, 3j k y k G 7 ifc } (1-5) 

with a measure dM 1 dM 2 such that if (x±, . . . , xn) are the eigenvalues of Mi and (yi, . . . , y^) are 
the eigenvalues of M 2 , and x k G 7« fe and y k G jj k : 

dM 1 dM 2 :=det(K ikJl e Bi kVH) ( JJfo - x ,) J ( J]^ - Vj ) j dCA 2 J] ^ J] d W (1-6) 

\i>j / \«>i / i=l i=l 

Our goal is to compute the Schwinger-Dyson equations, also called loop equations in the 
context of random matrices, of the following matrix integral: 

Z:= f e-^mm+VtW-MM] dMi dM2 

1.1 Notations 

Define the following polynomials, which vanish at all the hard edges: 

s(x) := H(x - Xj) = £ s r x r (1-8) 

j=l r=0 
j=l r=0 

The resolvent, (which is a formal series in its large x expansion), is defined by: 

And, up to a shift by the potential, it is more convenient to use Y(x): 

Y(x) :=V((x) -W(x) (1-11) 

Then, for technical intermediate calculations, define the following expectation values, which 
are formal series in their large x expansion: 



N \ x - Mi y - M 2 

(1-12) 



x - Mi y - M 2 - Mi 

(1-13) 



c 



N \ x-M x y - M 2 



B(x) := LU^L KMi ) 
y ! N \ x — Mi 



N \ x-Mx 



(1-14) 



;i-15) 



;i-16) 



;i-i7) 



K 2 3-1 

D(x) := xB(x) +J3E S - M ( TlM ^ J ) B >W (US) 

Notice that U (x, y) and U (x, y, x') are rational functions of y, and P(x, y), is a rational function 
of both x and y, whose poles are the same as the poles of V[ and V^. 



2 The loop equations 

Loop equations, i.e. Schwinger- Dyson equations, merely express the fact that an integral is invari- 
ant under a change of variable. Schwinger equations have been extensively studied in the context 
of field theory, and have been of great improtance for the study of matrix models j^j. For the 
polynomial 2-matrix model, loop equations were first announced to give algebraic equations by 
M. Staudacher jTH], and computed precisely for general potentials, including 1/N 2 corrections by 
[T2*| EE QUj- They have been solved recentely by [8. following the method of [Zj. Here, we 
follow the method of [T3|. 

2.1 Generalities 

2.1.1 Changes of variables 

Consider a change of variable in H 2 (T), of the form: 

M. - M[ M. I ef{M u M 2 ) + ef\M u M 2 ) + 0(|e| 2 ) (2-1) 



Remark 2.1 Under such a change of variable, the eigenvalues of M[ are no longer on paths ji h , but on 
some paths 7^ , which, for e small enough, are small deformations of ji k , and thus the integral is indeed 
unchanged. This is true only if we can deform the contours, i.e. not at the extremities. At poles, the 
integrand vanishes, so that the integral is still invariant. At hard edges, the integral is invariant only if 
M[ = Mi, i.e. f(M\,M2) vanishes at hard edges. 



To order 1 in lei, one has: 



exp 



N 

-Tr (Fi(M0 + V 2 (M 2 ) - M[M 2 ) 



(l-eS(f)-eS(f) + 0(\e\ 2 )) exp 



N 



Tr (Vi(Mi) + V 2 {M 2 ) - M ± M 2 ) 



with 



N 



S(f) = - Tr (VKM,) - M 2 )f(M 1 , M 2 ) 



The measure dMi is multiplied by a Jacobian J, which we expand to order 1 in |e|: 

dM[ = J(f) dM 1 = (1 + ej(f) + ej(f) + 0(e 2 ))dM l 
Loop equations are obtained by writing that the integral is unchanged, i.e.: 

Z := / dMidM 2 exp 



dM{dM 2 exp 



N . 



— -Tr (VI (Mi) + V 2 (M 2 ) - M 1 M 2 ) 
N 

--Tr (Vi(M() + V 2 (M 2 ) - M[M 2 ) 



dMi dM 2 (1 + e (J(f) - S(f)) + 1 (J(f) - S(f)) + 0(|e| 2 )) 



exp 



N 

--Tr (^(Mi) + V 2 (M 2 ) - MiM 2 ) 

v 



(2-5) 

Since this must be true for all argument of e, one must have: 
0= / dM 1 dM 2 (J(f)-S(f)) exp 



N 



— Tr (y x (Mi) + V 2 (M 2 ) - M X M 2 ) 

v 



i.e. 



(J(f)) = W)> 



(2-2) 



(2-3) 



(2-4) 



(2-6) 



;2-7) 



This equation is called a loop equation. 

Let us emphasize again that eq !2-5|. and thus eq !2-7l hold only if / vanishes at the hard edges. 



2.1.2 Split and merge rules 



The rules to compute J~(f) are called split and merge rules and can be found in the litterature, 
for instance in [Hj. Notice that J is linear, and if / is a product, J(f) can be computed by the 
chain rule. Thus, it is useful to determine J{f) for some particular /. For any two matrices A 
and B, one has: 



Split rule: 



or equivalentely: 



Merge rule: 



or equivalentely: 



J I A K^B] = TtA ^Tr K^B 



x - M l 



x - Mi x - Mi 



fc-i 



J (AM^B) = TrAM( TrMf'^'B 
3=0 



J ATi 



B 



Tr — - A l — B 



x - Mi J J x - M 1 x-Mi 

k-l 

J (A Tr (M x fe B)) =J2 TrM i A M^ 3 ' 1 B 



(2-8) 



>9) 



(2-10) 



(2-11) 



3=0 



2.2 Equations 

In this section, we write the change of - variable, i.e. function /, and the corresponding loop 
equation of type eq !2-71 We always write the contribution of the Jacobian J in the LHS and the 
contribution of the action S in the RHS. 
• computation of B{x): from 



f(M 1 ,M 2 ) 



s(M 1 
s(x)- 



s(Y(x)) - ~s{M 2 ) 



x - Mi Y(x) - M 2 
1 s(Y(x)) - s(M 2 ) s(x) - s(Mi) s(Y(x)) - s(M 2 ) 



x - Mi Y{x) - M 2 
(which indeed vanishes at all hard edges) we get: 



x-Mi 



Y{x) - M 2 



(2-12) 



s(x)W(x)±- (tt^— *~(^)) - W 
1 ' 1 >N\ x-M l Y{x)-M 2 

K J N 2 \ x-M 1 x-Mi Y(x) - M 2 



r =0 i=0 ^ 



r— 1— i 



Tr 



x i - M\ s(Y(x)) - s(M 2 ) 
x - Mi Y(x) - M 2 



- iv \ Tr ^Mi" (v; (Ml) - Ms) y(x) - m 2 / 

(2 - 13) 

i.e. 



= s(x)s(Y(x))W(x) 

I / V>{x) - V{(M X ) ~s(Y(x))-~s(M 2 ) \ 

1 J iV \ x-Mi -M 2 / 

* / m s(x) -s(Mi) ,,,, 



iV \ a; -Mi v 1V 7 ' - M 2 

5(y(x)) - 5(M 2 ) 
iV 2 ^ ^ r \ 1 x - Mi y(x) - M 2 

v ; iV 2 \ a; -Mi a; - M x y(ar) - M 2 / c 

(2 - 14) 
• computation of Bk(x): from 

/(Mi,M 2 ) = S (M0^— y( f\ fc ~ M2fc 
V i; x-Mi y(ar)-M 2 

1 F(x) fc - M 2 fe _ s(x) - s(Mi) y(x) fe - M 2 fc 
S ^x-Mi y(rr)-M 2 x - Mi y(rc) - M 2 

(which indeed vanishes at all hard edges) we get: 

s{x)W{x)— ( Tr 



N \ x-M 1 Y(x) - M 2 

t 2 / i i r(x) fc -M 2 fc \ 

+ s Wjp \ Tr ^M^ Tr ^M^ y^) - M 2 / c 

^ 4^ 4) \ x - Mi y(x) - M 2 

l/^-^MMO-M)^-^ 



(2 - 16) 



iV \ x-Mi 1V y ' y(x) -M 2 



i.e. 



B k (x) 



Y(x) k s{x)W{x) 



Tr 



V{{x) - ^(s) 



M 2 fc 



a; 



Mi 



y(x) - m 2 



t I ■ s(x)-s(M 1 ) Y(x)*-M* \ 
~N \ Tr x-M 1 {Vl (Ml) ~ M2) Y(x) -M 2 ) 

N\ x - Mi V 

N2 hh \ *- M ^ n*)-Mj 

i 2 / 1 1 Y{x) k - M\ \ 

~ s Wjp \ Tr ^M; Tr ^Mi y(x) - M 2 / c 

(2 - 17) 
• Computation of D(x): 



D(x) = xB(x) 

s-l 

Vs. 

AT 



s=0 j=0 

xs(x)5(F(x))Vr(a;) 

rqM l / T V i( x ) ~ v (( M i) §(Y(x))-s(M 2 ) \ 
XSlJ iV\ ir x-M l Y(x)-M 2 / 

s=o i=o \ 1 v 7 2 

-YYh- (Tr Mr 1 -) - / Tr S{x) - 3{Ml) {V[{ Ml ) - M 2 ) ^ = M ° 
N X 2 1 N\ x-Mi K U ' ' Y(x)-M 2 

s=0 i=0 r=0 i=0 \ i V ; z 

_ xs(x) il /-,>_!_ ggWLM\ 

V ; N 2 \ x-M 1 x-M 1 Y(x) - M 2 I c 

s=o j=o \ 1 1 v 7 z 



s=0 j=0 N 1 

(2 - 18) 

• computation of A(x): By doing a change of variable on M 2 of the form: 

f(M u M 2 ) = s ( Ml )—^s(M 2 ) 
(2-19) 



we have: 



2 ^2 s-1 



i.e. 



AT2 

- xB(x) + (Tr s(Mi) s(M 2 )> 

(2 - 20) 



= £>(> 



--(Tr S (M05(M 2 )) 

+4 ^ ? * / Tr -i*^- M| Tr M*' 1 ' 3 
N2 \ x - M 1 2 2 



s=0 j=0 

(2-21) 
Main computation: from 

1 s(Y(x))V 2 \Y(x)) - 5(M 2 )\/ 2 '(M 2 ) 



f(M 1 ,M 2 ) = s(M 1 ) 



= six 



x - Mi y(x) - M 2 

1 5(y(x))^'(^(^)) - 5(M 2 )y 2 '(M 2 ) 



x - Mi y(x) - M 2 

s (x) - s(Mi) 5(y(x))y 2 , (y(x)) - §(M 2 )y 2 '(M 2 ) 



we get: 



x - Mi Y(x) - M 2 



s(x) l W(x) I Tr 1 5(y(x))^(y(x))-5(M 2 )^(M 2 ) \ 
s WjV^ W \ x-Mi Y[x)-M 2 / 

+s(x)— /Tr 1 Tr 1 " \ 



(2-22) 



N 2 \ x-M 1 a: -Mi y(x) - M 2 / c 

iv^Z^Z^ 5 ^™! ir 7^M7 y(x) - m 2 / 

r=0 i=0 



(2 - 23) 

i.e. 



1 _J_ W <*> - M 2) 



* 2 / i i 5(y(x))^(y(x)) - g(M 2 )\/ 2 -(M 2 ) \ 

W 7V 2 \ x-Mi x-Mi y(x)-M 2 / ( 

s(x)s(Y(x))V 2 \Y(x))W(x) 

_ a ( x )l / Tr^'^ " Vl ' (Ml) W^W*)) " 5(M 2 )V 2 '(M 2 ) \ 



(2 - 24) 

i.e. 



iV \ x - Mi - M 2 / 
* L ^)-«W , 5(y(a;))^(y(x)) - 5(M 2 )V 2 '(M 2 ) * 
~iV \ Tr x-M x ( ^ l(Ml) ~ M2) Y{x) - M 2 , 

A 1 x-Mj y(x)-M 2 , 



^ / Tr 1 Tr 1 g(y(x))^(y(x))-g(M 2 )^(M 2 ) \ 

1 ' N 2 \ x-M x x-M x Y(x)-M 2 / £ 

+ 4 V V 5 S / Tr M 2 Tr M^^A 

iV 2 ^^ \ re -Mi 2 2 / 

s=0 j=0 x 1 1 c 

s(x)g(y(a;))y 2 / (y(x))^(a;) 



N \ x .1/, ~ ~ V(.r) - M 2 / 

f /^ s(x)- S (Mi) 5(y(x))V^(y(x)) - g(M 2 )\/ 2 '(M 2 ) * 

~N \ Tr x-Mi (Vl(Ml) " M2) y(x) - M 2 , 

+ilvy% /w- 1 - t / " M * !gMMg^j - gWjgW N 

^ 2 ^f^ r \ 1 x-Mi y(x)-M 2 , 

+l(Tr S (Mi)5(M 2 )) 
-D(ar) 

(2 - 25) 

Now, inserting the value of D(x), and after tedious but straightforward computations, we 

s(x)s(Y(x))— I Tr - Tr - M^MhlM^l \ 

St J4 i ,J JV 2 \ x-M 1 ir x-Mi y(x)-M 2 / c 



\ x-Mx y(x)-M 2 / c 

^"2 s-1 2 / I \ 

^mEE^™^— m^) 

s=0 j=0 x 1 / c 

s=0 j=0 r=0 i=0 \ iv/ 

s(x)s(Y(x)) m(Y(x)) - x)(Vl(x) - Y(x)) - P(x, y) + t) 



r=0 i=0 

_t_ / s(x) - s(M 1 



(2 - 26) 



N\ x - Mi 

(\/ 1 / (M 1 )V 2 '(M 2 ) - il^'W - MxVYCMx) + MxM 2 ) S( ^j " ^ 

g / s(x) - s(M 1 ) ^ g(y(x)) - g(M 2 ) \ 
iV 2 \ x - Mi - M 2 / 

4 E X> <™-> <i£$ «<*> - mo !03g>z3*> ) 

r=0 i=0 N 1 V 7 Z ' 

4EEEE- ^> (^^g 

r=0 i=0 s=0 j=0 \ i \ J z 



which we write as follows: 



s(x)s(y) E(x,Y(x)) = j^L(x) 



(2-27) 



where E(x,y) is a rational fraction of both x and y with poles at the poles of V{ and V 2 , and at 
the hard edges: 



E(x,y) 

(Vl(y)-x)(V((x)-y)-P{x,y)+t 
V{{x)-V[{Mi) 



1 t 



s(y) N 

K 2 s-1 



Tr- 



x - Mi 



y-M 2 

N 

s=0 j=0 

1 t I six) - s(Mi) /Tr//W s 

Tr (V^ (Mi) - M 2 ) 



six) N 



x - Mi 



x - Mi j/ - M 2 



l t 



2 Ki r ~ l I 



\ 7 r=0 i=0 x 



TrM[- 1 - i Tr 



y-M 2 
x'-Mj V^y)-V 2 '(M 2 ) \ 
x - Mi y - M 2 / 



I t s(x) - S (Mi) 



s(x)%) iV \ x - Mi 

(^'(MOl/^M,) - M 2 V 2 \M 2 ) - MxF^Mi) + .\ /,.!/,; 

1 - g(Mi) ^ %)-g(M 2 ) \ 

s(x)s(y) N 2 \ T x - Mi ' y - M 2 / 



%)-g(M 2 ) \ 
y - M a / 



1 t 2 Kl r ~ 1 



s(x)s(y) N 2 r ( . ,_ () 

1 , 2 ^2 S-1 



x -"M" ( ^ (Ma) " Ml) y-M 2 / 



s(x)s(y) N 2 s _ {] ( . 

3 r— 1 s— 1 



.s, (Trl/r'-Q (Tr^4^ (^(MQ - M 2 ^ ^ " M * 



x - Mi 



1 



t 6 



s(x)s(y) N 3 



EEEE 5 ^ (TrMr 1 -) (TrMr 1 -) Tr 



r=0 i=0 s=0 j=0 



y-M 2 

x l - M\ y j - Mi 
x - Mi y-M 2 



(2 - 28) 

and in the right hand side of eq !2-27| we have: 



Tr 



1 



1 



, ( r . vg(y(x)) - ^(m 2 ) \ 

x - Mi x - Mi F(x) - M 2 / 



Tr 



1 



Tr- 



s(Y(x)) \ x- Mi x- Mi 



(V'(M) M .) ~ S(X{X)) -~ S{M2) \ 
{V 2 {M 2 )-Mi) Y{x) _ M2 ) 



L_ / Tr i g(y(x)) - g(M 8 ) \ 

s(y(ar)) \ x - M x - M 2 / c 

1 ^ 2 5-1 / 1 



K 2 s-1 



t / „ i j\ / 1 m 1 y(a;V'-M| 

/TV. fl/f s ~ x ~ J\ / TV. TV. V / ^ 



* y y ~ Ss 1 ( Tr MT^') ( Tr^^- Tr 



s(Y(x)) N x 1 1 \ x - Mi x - Mi - M 2 

iCi r-l K2 s-1 

s(s)s(y(a;)) 



+ 77^Wwvri S SrSt 



r=0 i=0 s=0 j=0 

1 ( Tr mth) / Tl . Mr i- Tr £^M MiM ' 

K, r-l 



(2 - 29) 



2.3 Examples 

Equation eq !2-27l looks rather terrible, but it is actualy very simple to use. Let us illustrate it on 
simple examples. 

2.3.1 No hard edges 

If there is no hard edges, we have s(x) = 1 and s(y) = 1, and the loop equation becomes: 

(TO -z)(y{(x)-Y(x))-P(x, r W)+ «4 Tr-^- ™I^ ) c 

(2-30) 

which is the well known loop equation of the 2-matrix model with polynomial potentials [T%] . 

2.3.2 1-matrix model 

Consider the 1-matrix model with weight e~^ TTV ^ M \ It is equivalent to a 2-matrix model where 
M 2 is gaussian, i.e. V 2 (?/) = y, s(y) = 1, V'{x) = V{(x) — x. The equation eq !2-27l becomes: 

\ , WT ,„ , . t / m V'fo) - l"(M)\ 
= (y-x)(y'(x) +x-y) - — (Tr— ^ K -^\ 



N \ x-M j 



1 t / Ti s(x)-s(M) v , (M) 



s(x) N \ x-M 

'y 7 / . . - }[' 



2 ai r-i , 

i-^-VV Sr TrM r_1_i Tr: 
s(x) N 2 ^ ^ \ x-M 

\ 7 r=0 i=0 x 



(2-31) 

and in the right hand side of eq 12-271 we have: 

Since E(x,y) is quadratic in y, this equation defines an hyperelliptical curve. 

2.3.3 1-matrix model with only one hard edge 

In particular, consider s(x) = (x — a), we have: 

E(x, y) = (y- x)(V'(x) +x -y)-±(rr ^^I^ ) - < Tr V '( M )) ( 2 " 33 ) 

2.3.4 1-matrix model with only two hard edges 

In particular, consider s(x) = (x — a) (x — b), we have: 



E(x,y) = (y — x)(V'(x) + x — y) — — ( Tr 



t V'jx) - V'{M) \ 



N\ x-M J 



1 ' (Tr (x + M - a - b) V\M)) ' ' 



(x — a)(x — b) N ' (x — a)(x — b) 

(2 - 34) 

2.3.5 2-matrix model with only one hard edge 

Consider s(x) — (x — a) and s(y) = 1, we have: 

E(x,y) = (Vi(y)-x)(Vl(x)-y)-P(x,y) + t 

1 ( /tv W( mo-m 2 ^-™ 



x-aN \ y iV ' ' y-M 2 
(2 - 35) 



and in the right hand side of eq 12-271 we have: 



L(x) _ (Tr ^ Tr 1 ^M)-v?(M 2 ) 



x - Mi x-Mi Y(x) - M 2 
(2 - 36) 



2.3.6 2-matrix model with only two hard edges 

Consider s(x) = (x — a){x — b) and s(y) = 1, we have: 

E(x, y) = (V&y) - x)(V[(x) -y)- P(x, y) + t 



1 ' /iK. + Jt-.-mw-jft)^-™ 



s(x) N \ v ,uv ' 1 y-M 2 

I 1 / V>(y)-V'(M 2 ) \ 
s(x)N\ y-M 2 J 

(2 - 37) 

and in the right hand side of eq !2-27[ we have: 

L(x) = / T r 1 Tr 1 ^ ~ WW \ 
K) \ x-M x x-Mi F(x)-M 2 / c 

(2 - 38) 

2.3.7 2-matrix model with only one hard edge in x and one hard edge in y 

Consider s(x) = (x — a) and s(y) = (y — b), we have: 

E(x,y) = (V 2 (y)-x)(Vl(x)-y)-P(x,y)+t 



s(y) iV \ x - Mi 

1 (Tr(Vl(M 1 )V 2 \M 2 )-M 2 V 2 \M 2 )-M 1 V;(M 1 ) + M 1 M 2 )) 



s(x)s(y) N 
t 2 _ 

s(x)s(y) 



(2 - 39) 

and in the right hand side of eq 12-271 we have: 

L w = \ Tr - — rr^ 



x - Mi x - Mi F (x) - M 2 

(2 - 40) 



s(r (x)) \ x — Mi x — Mi 



3 Large AT limit, algebraic curve 

In the large N limit, eq 12-271 reduces to an algebraic equation: 

E(x,Y(x)) = 



(3-1) 



One should notice that in the large N limit, there is the factorization jB]: < TrTr >=< Tr >< 
Tr >, so that x and y play symmetric roles [T7j . 

We see that E(x, y) has poles only at the poles of V[ and zeroes of s in x and at the poles of 
V% and the zeroes of s in y. 



3.1 Behaviour near hard edges 

near an hard edge x — > X iy such that s(Xi) = 0, we have: 



Thus: 



six) iV x 1 / N \ x - Mi 

v ' r=2 1=1 N 1 

+finite (3-2) 



1 



ResF (x)dx = — ; — - — ( Tr 
(3-3) 

Hard edges are at the same time poles of Y(x), and zeroes of dx, so that Y(x)dx is finite. 

3.2 Behaviour near poles of the potential 

Near a finite pole £ of V{(x), we have: 

y(x) ~ V((x) - W{£) + 0(x - (3-4) 



thus 



KesY(x)dx — KesV((x)dx (3-5) 



Near a pole at 00, we have: 

t 



thus 



Y(x) ~ y/(x) - - + 0(1/ x 2 ) (3-6) 



ResY(x)dx = t (3-7) 



3.3 Determination of the algebraic equation E(x, y) 

So far, we know that E(x, y) is a rational function of x and y, we know its form, and its behaviour 
near poles, but most of the coefficients are not determined by the loop equations. 

The remaining coefficients of E(x,y), as usual, are determined by extra requirements, which 
depend on how the matrix model is defined, i.e. on the purpose for which we introduce the matrix 
model. 

The two most frequent definitions of the matrix model are (we mainly follow the presentation 



3.3.1 Case of the convergent matrix model 

In this case, the matrix model is defined by the convergent integral eq 11-71 For generic potentials 
and hard edges, -^InZ has a large N limit F Q , but has no 1/N series expansion. The resolvent, 
and thus the function Y(x) also, has a large N limit, but no 1/N expansion. This fact can be 
understood from the work of 0j. 

The large N limit of Y(x) obeys an algebraic equation E(x, Y(x)) = 0. 

Now, consider an arbitrary E(x, y) satisfying the correct behaviours near poles. It gives a 
function Y(x), and from it one can compute the free energy Fq (see the formula in |2j). The Y(x) 
which is the large N limit of the resolvent is the one for which dt(Fo) is minimal (see [THll^]). That 
implies in particular that for any contour C on the algebraic cmveE(x, y) = 0, one has: 



On an algebraic curve of the type E(x, y) = 0, there are at most 2 x genus independent irreducible 
cycles, (plus contours around poles), and one can check that the number of constraints of type 
eq 13-81 exactly matches the number of coefficients of E(x, y) not determined by the pole behaviours. 
Thus, condition eq 13-81 determines E(x,y) completely. In case there are several solutions, one 
determines a unique one, by choosing the absolute minimum of F Q . 

3.3.2 Case of the formal matrix model 

The formal matrix model can be defined in a combinatoric way, as a formal series, generating 
discrete surfaces. The formal expansion is obtained by expanding the matrix integral eq ll-71 with 
the Feynman graph technics EE E] • 

In that model, InZ has a 1/N expansion by definition, as well as the resolvent, and all ex- 
pectation values. The resolvent of the formal matrix model, is thus obtained by solving the loop 
equation eq !2-27| order by order in 1/N 2 . 



of [El): 




(3-8) 



That model, in addition to the potentials and hard edges, depends on a "vacuum" around 
which the Feynman expansion is performed. This vacuum is characterized by a set of "filling 
fractions", as follows (see (THj for more details): The potential V\(x) + V^y) — xy has a certain 
number of extrema, which are given by the algebraix equation: 

V((x)=y, Vi{y) = x (3-9) 

i.e. 

V&V&x)) = x (3-10) 
Let us call K the degree of that algebraic equation, i.e. the number of its solutions: 

(x l ,y 1 ),...,(x K ,y K ) (3-11) 

The eigenvalues of matrices M±, Mi which extremize Tr(Vi(Mi) + V2(M2) — M\M<i) must be among 
the K solutions described above, or can be trapped on contours stopping at hard edges. The filling 
fractions (ei, . . . , tK+K 1 +K 2 ) represent the number of eigenvalues equal to each solutions of eq !3-101 
i.e. there are iVei eigenvalues of M\ equal to x\ y .. One must have: 

K+K^+K 2 

E e * = 1 ( 3 - 12 ) 

The average number of eigenvalues of M\ in the vicinity of a point (x, y), is a contour integral of 
the resolvent, along a contour which surrounds {x,y): 

2mNei = - i W{x)dx = i Y(x)dx (3-13) 
J A J Ai 

and it can be non-zero only on irreducible cycles of the algebraic curve. 

Thus, in the Formal matrix model, a set of irreducible cycle contour integrals are fixed. One 
can verify that the number of filling fractions matches exactely the number of coefficients of E(x, y) 
not determined by the pole behaviours, Thus, eq !3-13l is sufficient to determine completely the 
rational function E(x,y). 

4 Conclusion 

The purpose of this article was to write down the loop-equations for the so-called semiclassical 
2-matrix model. We find, that the loop equation becomes an algebraic equation E(x,y) = in 
the large iV limit, with poles at the poles of the potentials, and at the hard edges. The hard edges 
are such that the resolvent has a simple pole, and dx has a zero, so that the differential form ydx 
is regular. 



The loop equation determines the form of the algebraic equation E(x,y), but does not de- 
termine its coefficients. The coefficients are determined by extra asumptions, related to which 
definition of the matrix model is consiered. In the formal matrix model, the A cycle integrals of 
ydx are fixed parameters of the model, and that determines E(x,y) completely. In the conver- 
gent matrix model, the real parts of both A and B cycle integrals of ydx must vanish, and that 
determines E(x,y) completely. 

Let us also remark that the function eq \'2-27\ E(x. v) is unchanged under the exchange x <-> y, 
V\ <->• V2, s <-» s, which is the generalization of Matytsin's duality property [T7j : 

X(Y(x)) = x (4-1) 

The consequences of that algebraic equation, can then be studied. This is done for instance 
in One can also expect to generalize the works of jOj, or or |H], i.e. the computation of all 
correlation functions and their 1/N 2 expansion, and further, compute the expansion of the free 
energy jTUHH]. 
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